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Abstract 

It has been experimentally observed that weakly conducting suspended films 
of smectic liquid crystals undergo electroconvection when subjected to a large 
enough potential difference. The resulting counter-rotating vortices form a 
very simple convection pattern and exhibit a variety of interesting nonlin- 
ear effects. The linear stability problem for this system has recently been 
solved. The convection mechanism, which involves charge separation at the 
free surfaces of the film, is applicable to any sufficiently two-dimensional fluid. 
In this paper, we derive an amplitude equation which describes the weakly 
nonlinear regime, by starting from the basic electrohydrodynamic equations. 
This regime has been the subject of several recent experimental studies. The 
lowest order amplitude equation we derive is of the Ginzburg-Landau form, 
and describes a forward bifurcation as is observed experimentally. The coeffi- 
cients of the amplitude equation are calculated and compared with the values 
independently deduced from the linear stability calculation. 
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I. INTRODUCTION 



Although spatio-temporal pattern formation is ubiquitous in nature, there are relatively 
few systems which are amenable to both accurate experimental study and first-principles 
weakly nonlinear analysis. [|l[] The classic examples involving fluid mechanical instabilities 
are Rayleigh-Benard convection and Taylor- Couette flow. The results of perturbation theory 
based on the Navier-Stokes and heat equations are in good agreement with high precision 
experiments in the weakly nonlinear regime of these two instabilities. A more complex 
example is electroconvection in nematic liquid crystals due to the Carr-Helfrich mechanism. 
H Here, remarkably good agreement has been achieved in spite of the complexity of the 
problem. However, in other cases either the materials cannot be sufficiently characterized 
or the underlying physical equations are not understood well enough to allow quantitative 
comparisons between observations and theory. 

Electroconvection in suspended smectic films is a promising new experimental system for 
quantitative studies of spatio-temporal pattern formation. |4]-|7j When a thin, suspended film 
of smectic liquid crystal is subjected to a sufficient potential difference, a charge separation 
arises which drives the film into convection. The flow pattern just above onset is sustained 
by the electric field acting on charges which develop near the free surfaces of the film. These 
charges are simply a consequence of the electrostatic boundary conditions which must be 
satisfied by the fields inside and outside of the film. || Fig. [I] shows a schematic of the 
experimental arrangement. This source of charge is distinct from that due to the Carr- 
Helfrich mechanism which drives bulk electroconvection in certain nematics. || In that case, 
the charge generation mechanism involves an essential coupling to the director orientation. In 
experiments on smectic A liquid crystal films in which the director was perpendicular 

to the film, no orientational effects were observed, indicating that the flow remained isotropic 
in the film plane. Recent experiments on smectic C films showed convection and flow 
alignment of the projection of the director in the plane of the film, but were not consistent 
with the Carr-Helfrich mechanism. These were likely driven by the mechanism discussed 
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here, with the flow alignment a secondary effect. It has, however, been hypothesized flO 



that the two mechanisms might coexist in some smectic C materials. The mechanism we 
discuss here is presumably also responsible for convection observed in thin, suspended films 



of isotropic fluids and in nematics in certain regimes JTTJ. These cases involve substantial 
three-dimensional effects, however, because they lack the smectic layering which has the 
effect of restricting the flow to the film plane. In what follows, we consider only very thin 
isotropic films, relevant to the case of smectic A, on which most of the experiments have 
been performed. f§-[7| 

A theoretical model of the onset of electroconvection in suspended films was introduced 
by Daya, Morris, and de Bruyn. || The film was represented as a weakly conducting, two- 
dimensional, isotropic fluid. To find the electric fields and charge densities which drive 
convection, the electrostatic potential was determined. The three-dimensional electrostatic 
equations effectively constitute a nonlocal coupling between the in-plane fields and charge 
densities which appear in the two-dimensional Navier-Stokes and charge continuity equa- 
tions. This extra coupling formally distinguishes the resulting equations from those of ther- 
mal convection in the Rayleigh-Benard problem, although some interesting similarities re- 
main. The value of the critical wavenumber from the linear stability analysis || is in good 
agreement with experiment HH|7]]. 

The purpose of this paper is to present a weakly nonlinear analysis of the electrohydro- 
dynamic equations given in Ref. |J. The multiple-scales perturbation theory employed in 
our treatment is similar to that given in Ref. |]] for Rayleigh-Benard convection. Although 
there are important physical differences between these two pattern forming instabilities, the 
resulting amplitude equation for both problems is of the Ginzburg-Landau form 

r d t A = eA + e d 2 x A-g A\A\\ (1.1) 

where e is the control parameter, which depends on the applied electric potential, and 
A(x,t) is a slowly varying amplitude. The coefficients r , £ , and g are compared with 
those previously obtained by other methods. In particular, r and £o are found to be in 



good agreement with the values determined from the linear stability analysis of Daya et al. 
H Mao, de Bruyn, and Morris have experimentally measured all three coefficients. |7| The 
experimental value of £o is in reasonable agreement with theory. Quantitative comparison 
of ro and go with theory is difficult at the present time due to large uncertainties in the 
conductivity and viscosity of the liquid crystal, which are required to non-dimensionalize 
the experimental results. 

Determining the amplitude equation constitutes a first step towards understanding the 
fully nonlinear regime beyond the onset of electroconvection. For small wavenumbers near 
threshold, the stability of solutions of Eqn. determines the regions of control parameter- 
wavenumber space where the vortex pattern itself is stable. |IJ For example, one expects 
such a one-dimensional pattern to exhibit a long-wavelength instability due to the Eckhaus 
mechanism which restricts the range of stable wavenumbers available to the pattern. 

The amplitude equation can also be used to study how the ends of a finite-length film 
affect the range of stable wavenumbers. This wavenumber selection mechanism was first 
investigated by Cross, Daniels, Hohenberg, and Siggia [113 in the context of Rayleigh-Benard 
convection in finite containers. End-selection was observed experimentally by Mao, Morris, 
Daya, and de Bruyn |)J in electroconvection patterns in smectic A films. It is possible to 
extend the present theory to determine the Eckhaus and end-selection stability boundaries 



14| , but this is beyond the scope of the present paper. 

This paper is organized as follows. The linear stability analysis of the electroconvective 
instability is briefly reviewed in Section ||. The amplitude equation is determined in Sec- 
tion |TJ. Section [IV] compares the results of this theoretical investigation with previously 
obtained observations, and contains a brief conclusion. 



II. LINEAR STABILITY ANALYSIS 

In this section, the physical model describing electroconvection in a thin film is presented. 
The linear stability analysis of the relevant equations is concisely reviewed. Further details 
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are given in Ref . || . Note that we change some of the notation of Ref. || to simplify the 
presentation of this paper. 

The film is treated as a two-dimensional conducting fluid in the xy-pl&ne, with areal 
material parameters p s = sp, r] s = sr), and a s = sa, where s, p, r], and a are the film 
thickness, bulk density, bulk molecular viscosity, and bulk conductivity, respectively. The 
coordinate system is shown in Fig. |l|. The film is assumed to be infinite in the x-direction 
and bounded between — d/2 and d/2 in the ^/-direction. We only consider the thin film limit 
where s/d — > 0. Two electrode configurations are analyzed. In the 'plate' geometry, the film 
is suspended between two thin sheet electrodes which fill the rest of the xy plane, whereas in 
the 'wire' geometry, the film is suspended between two thin line electrodes which are along 
the x direction. In both cases, the electrode at y = —d/2 is fixed at a potential of — V/2, 
and the electrode at y = d/2 is at a potential of V/2. 

The Navier-Stokes equation 



Ps 



^ + (u-V.)u 



-V S P S + r] s V s 2 u + qE s , (2.1) 



describes the fluid flow of the liquid crystal, where V s = (d x ,d y ,0), P s (x,y), q(x,y), and 
E s (x, y) are the two-dimensional gradient operator, two-dimensional pressure field, surface 
charge density, and electric field in the film plane, respectively. The incompressibility of the 
fluid implies that 

d x u + d y v = , (2.2) 

where u(x,y) and v(x,y) are the x and y components of the two-dimensional velocity field 
u. The pressure field is eliminated from Eqn. |2.1| by applying the curl operator. Taking the 
curl of Eqn. BTT] twice, using Eqn. |2]2|, and selecting the y component gives: 



- p s d t V 2 s v + p s d y (V s ■ [(u ■ V s )u]) - p s V 2 s [(u ■ V>] = 
-r) s Vtv + (d 2 x q)(d y *\ z=0 ) + {d x q)(d xy *\ z=0 ) ~ (d xy q)(d x *\ z=0 ) - (d y q)(d 2 x V\ z=0 ) , (2.3) 

where the electric potential \l/(x, y,z) is related to the in-plane electric field via E s (x, y) = 
— V s ^(x, y, z)\ z=0 . The three-dimensional Laplace equation, 



V 2 ^ = 0, 



(2.4) 



specifies \I/ in the half space z > with appropriate boundary conditions in the :q/-plane, 
where V = (d x ,d y ,d z ). The surface charge density depends on the discontinuity in the 
z— derivative of \I/ across the two surfaces of the film: 



q = -eo 

= -2e 



dz 
dz 



=0+ 



oz 



2=0- 



(2.5) 



=0+ 



where eo is the permittivity of free space. 

The motion of charge is governed by the charge continuity equation 

dq 

■jfi = -V s - (qu + a s E s ), 



(2.6) 



which includes contributions from both the conductive a s Ei s (x, y) and convective 
q(x, y)u(x, y) current densities. Diffusion of charge in the plane of the film can be neglected. 

Equations |2.3|42.6| constitute the basic electrohydrodynamic equations; the electrode ge- 
ometry enters into the boundary conditions on The solution of these equations in the 
"base state" below the onset of convection has u^(x,y) = 0, with q^°'(y) and ty(°'(y, z) 
satisfying the electrostatic boundary value problem given by Eqns. |2.4| and |2.5| . To examine 
the stability of the base state, we introduce the perturbed quantities 



(i> 



u = + u 



E s = Ef' + EW 



(2.7) 
(2.8) 
(2.9) 



where = E y ^y and E^(x,y) = E x ^\x, y)~k + E y ( - 1 \x,y)y. Quantities which have 
dimensions of length, time, charge density, and electric potential are nondimensionalized 
by d, €od/a s , eoV/d, and V, respectively. Substituting the perturbed field variables into 
Eqns. |2~^-|2.6|, nondimensionalizing, and suppressing the superscripts, yields 
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- K%q + KQd 2 x S>\ z=0 = Kd x [(d x q)(d y S>\ z=0 ) - (d y q)(d x V\ z=0 )} 
+V~ 1 {d t (V 2 s v) - d y (V s • [(u • V>]) + V 2 [(u ■ V>]}, (2-10) 

-Qv + V^| 2=0 = d t q + u(d x q) + v(d y q), (2.11) 

!+(9»*)U - (9**)|z=o- = 0, (2.12) 

V 2 * = . (2.13) 

The dimensionless parameters TZ and V are analogous to the Rayleigh and Prandtl num- 
bers. We will henceforth consider only the limit V — > oo, as this is the case most relevant 
to experiments on real smectic materials 0, for which V ~ 10 — 100. The non-constant 
coefficient depends on the electrode configuration and is given by Q(y) = d y q(°\y), 

where q^(y) is the base state charge density. The variables v, q, and \I/ above represent the 
dimensionless perturbed functions q^\ and ^f^\ respectively, and satisfy the following 
boundary conditions: 

v(x,y = ±i) = (d y v)(x,y = ±~) = , (2.14) 



*(i,y = ±^0) = 0, (2.15) 



#(ar,j/,z) -> ,z -> ±oo . (2.16) 

In the plate electrode configuration, Dirichlet boundary conditions are employed on the xy 
plane, with 

q(x,y,0)=0, \y\>\- (2.17) 

In the wire electrode geometry, mixed boundary conditions apply such that 

d z V(x,y,z)\ z=0+ =Q, \y\>\- (2.18) 

In both cases the potential *$?(x, y,0) = ty\ z=0 (x,y), \y\ < 0, is specified on the film. 
Equations g7TO| - gTT3| can be expressed as 



CC = B 



(2.19) 



where 



/ 



Vf -Kdl TZQdl 
-Q V 2 











c 



o d z {.) 


v(x,y) 
q(x,y) 
y(x,y,z)\ z=0 
ty(x,y,z) 






=o+ - dz(-)\ z =o- 
V 2 



and 



B 



nd x [(d x q)(d y V\ z=0 ) - (d y q)(d x V\ z=0 )] 
d t q + u(d x q) + v(d y q) 



The linear stability problem is defined by 

CC = . 



V 



^2.20) 



(2.21) 



(2.22) 



(2.23) 



The neutral stability curve 1Z = 1Z c (k) is determined by substituting the normal mode 
solution 

f ~ ( ^ 
VkKV) 



C 



V 



e lKX = C K e l 



[2.24) 



q*(v) 

t«(j/,0) 

into Eqn. |2.23| . The following alterations have been made to the notation of Ref. |J: A(y) — > 
0(y) -> -> *«(y,0), and f%,2:) -> 2;). The variables u«, g K , and 

ty K are expanded as 



/ 
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oo 

(y) , (2-25) 

m=l 

oo 

QM = A m qKTn(y) , (2.26) 

m=l 

oo 

* K (y,z)= J2 A m* K m(y,z) , (2.27) 



m=l 



where v Km (y), q K m(y), and ^f Krn (y, z) satisfy the boundary conditions Eqns. |2.14| to |2TTE| . The 



linear problem is solved numerically in Ref. |§ by substituting C m , the even Chandrasekhar 
function fll6"l , for v Km (y) and then finding self-consistent solutions for q Km (y) and ty Km (y, z). 



III. DERIVATION OF THE AMPLITUDE EQUATION 

The multiple-scales approach is used to obtain the amplitude equation, which describes 
the slow temporal and spatial variation of the field variables. W^M The slow scales X = 
e l l 2 x and T = et are treated as independent of the fast scales x and t. We choose e = 
(TZ — 1Z c o)/1Z c o, where lZ c0 is the critical value of 71 at the minimum of the neutral stability 
curve TZ = 1Z c (k). The nonlinear equation describing the system, Eqn. |2.19| , is expanded in 
powers of e 1//2 as follows: 

C = C + e 1/2 C 1 + e£ 2 + ... , (3.1) 
C = t 1 ' 2 C + tC l +^ 2 C 2 + ... } (3.2) 
B = e 1 / 2 B + eB 1 +e^ 2 B 2 + ... , (3.3) 

and 

v(x, y) = e 1/2 v (x, y) + eu x {x, y) + e 3/2 v 2 (x, y) + ... , (3.4) 
q(x, y) = e 1/2 q (x, y) + eq x (x, y) + e 3/2 q 2 (x, y) + ... , (3.5) 
y, z) = e^oix, y, z) + e^ x {x, y, z) + e 3 / 2 ^ 2 (x, y,z) + .... (3.6) 

The differentials <9 X . and d t in the original equations transform as d x — > d x + e l ^ 2 dx and 
d t -> d t + ed T . 
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H/Cjll. y UeCOIIieS 




£oQ) — Bq , 


(3.7) 


£0^1 + CiCo = B\ , 


(3.8) 




(3.9) 



respectively. 

The most general solution of Eqn. 3/7, at order e 1 / 2 , is 



Co 



vo(x,y) 
Qo(x,y) 
-$ (x,y,z)\ z=0 



A (X,T) 



' -n N 
<?o(y) 

*oM) 
*o(^,z) 



o lK X 



+ c.c. 



(3.10) 



where Ao(X, T) is the amplitude function, Ko is the critical wavenumber which minimizes 
the function TZ c (k), and c.c. denotes complex conjugation. The functions vo(y) = v K (y)\ K=Ko , 
Qo{y) = 9k(i/)U=«o> and *o(2/ , z) = ^kG/, z)U=k are solutions of the linear stability problem, 
Eqns. [pi, |2~24l|2~27 . 

To solve the order e equation, Eqn. [3.8| , the relation 

^ — f \ iter ^ 



d_ 

dn 



C 



v K (y)e lK3 
q K (y)e iKa 

*«(?/,-2)e i 



/ 







(3.11) 



is employed. This is used to transform Eqn. 3.8 to 



£ Ci = Bi 



(3.12) 



which is solved by the ansatz 



Ci = 



qi(x,y) 

$i(x,y,0) 
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1 ^ 

vo(y) 


> 




1 \ 


Ql(y) 






A ie iK0X + c.c. 


> + 


9a(y) 












*S(y,o) 






K ^o(y, z)j 









|A,| 2 . (3.13) 



The variable Ai (X, T) is a second amplitude function. Substitution of Eqn. |3.13j into 
Eqn. [3.12| gives the following sets of partial differential equations 



(# - (2 Ko ) 2 )X + (2K yn c0 qi - (2K yn c0 Q^i 



2=0 



-2Kln c0 [q (d y ^o\ 



2 = J 



5i 



-Qui + (dy - (2k: ) 2 )^iU=o = iKoUoffo + v (d y q ) 

(d z n)\z=o+-(d z n)\ z=0 - = o, 

52 , o2 /r. 



(0„gb)*o|z=o] , (3-14) 
(3.15) 
(3.16) 
(3.17) 



and 



^ = , 

-(3^2 + 9^*2 1 *=o = -«'«o«o9o + iK ou* q + v (d y (fc) + Vo(d y q Q ) 



q e 2 + (d z %)\ z=0+ -(d z %)\ z=0 - = 
{dl + dl)n = 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



where the superscript * denotes complex conjugation. A vector C\ which solves Eqn. p. 12 
presented in Appendix A. The general solution at order e is 



is 



C 1 = £ - {(2 Ko )- 1 (2d x d x A )e^ x C' Q + c.c.} , 
where the prime denotes d K , and C' is given by 



C 







(d K v K (y))\ K=K0 


Qo(y) 




(d K q K (y))\ K=K0 


*o(y.«) 




(d K V K (y,0))\ K=Ko 






v (d K ^ K (y,z))\ K=K0 ) 



(3.22) 



(3.23) 
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The order e 3 / 2 equation, Eqn. [3.9| , can similarly be transformed to 



£0^2 



(3.24) 



by using Eqn. |3.11| . To establish the condition for the existence of a solution C2 of Eqn. |3.24 
the inner product 



(C b ,C a ) = (2tt/ko) 



2tt/kq /x 

- 1 1 dx 



dy {vl(x,y)v a (x,y) + q* b (x,y)q a (x,y) + %\ z=0 ty a \ z=Q } 



r27T / kq roc roc 

+ (2n/K y 1 dx dy dz {%(x,y, z)^ a (x,y, z)} , 

Jo J— 00 J— 00 



(3.25) 



is employed, where Q (i = a, b) is 



C; 



\ 



Vi(x,y) 
^i(x,y,z)\ z=0 

The adjoint operator £q is determined by integrating (C b , CoC a ) by parts: 

(C b , C C a ) = (£cA, C a ) + boundary terms. 



V 



(3.26) 



/ 



(3.27) 



The homogeneous equation CijZ a = with homogeneous boundary conditions, Eqns. |2.1 1 - 



2.18, will be referred to as the homogeneous "direct" problem. For this case, the left- 



hand side of Eqn. |3.27| is zero. Some of the boundary terms on the right-hand side of 
Eqn. |3.27 vanish because of the boundary conditions on the direct variables C a . By defining 
homogeneous boundary conditions for the adjoint quantities C b , the remaining boundary 
terms can be set to zero. This implies that 

C\C b = , (3.28) 



or more explicitly 



n c0 Qd 2 x v b + V 2 s q b + (d z y b )l 



V*u 6 - Qq b = 
-TZ c0 d 2 x v b + %\ z=0 = 
=o+-(3z*6)U=o- = 



V 2 ^ = . 



(3.29) 
(3.30) 
(3.31) 
(3.32) 
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Since Co 7^ mj, C is not self-adjoint. The corresponding operator in the Rayleigh-Benard 
problem can be made self-adjoint by an appropriate choice of the inner product. [[J This 
does not appear to be possible in the present case. The boundary conditions on the adjoint 
quantities are 



v b (x,y = ±i) = {d y v b )(x,y = ±~) = 



(3.33) 



%(x, y = ±i 0) = q b (x, y = ±-) = 



(3.34) 



-> ,2 -> ±oo 



(3.35) 



As in the direct problem, Dirichlet boundary conditions are employed on the plane in 
the plate electrode configuration with 

* b (x,y,0) = 0, M>~ (3.36) 

and mixed boundary conditions are applied in the wire electrode geometry with 

d z ^ b {x,y,z)\ z=0+ = 0, \y\>l- (3-37) 

In both cases the relation ^/ b (x, y, 0) = ^ b \ z =o(x, y), \y\ < 0, is specified on the film. 

The adjoint problem defined by Eqn. |3.28| with boundary conditions Eqns. |3.33| - |3.37| is 
satisfied by 



C h 



V 



-,IKQX n IKQX 



(3.38) 



VbKoW 

QbK (y) 

Equations |3.28| and |3.38| are analogous to Eqns. |2.23| and |2.24| in the linear stability problem. 
The details of the adjoint solution, Eqn. |3.38| , are given in Appendix B. 

Expanding the product (C b , CoC a ) for the inhomogeneous direct equation CoC a = Q, with 
homogeneous boundary conditions Eqns. [2.14H2.18| , leads to 



/ 
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(Cb, Q) = {C-\C b , C a ) + boundary terms. 



(3.39) 



The vector C b is specified by Eqn. |3.38| . The right-hand side of Eqn. |3.39| is zero due to 
Eqn. |3.28| and the boundary conditions on the direct and adjoint quantities. Equation |3.39 
therefore reduces to 



{C b , Q) = 



(3.40) 



where 



G 2 



+ C.C. + ... . 



(3.41) 



According to the Fredholm Theorem, Eqn. |3.40| is a necessary and sufficient condition for the 
existence of a solution of Eqn. |3.24| . [ |T5|j The remaining terms in Eqn. |3.41| do not contribute 
to Eqn. |3.40| since they do not depend on e ±tK ° x . Substituting Eqn. |3.38| and Eqn. |3.41| into 
Eqn. ^4D| yields 

(C w e iK0X , {G^e** * + c.c.) = ( (v b0 , g b0 , *&oU=o, ^w)e lK ° x , {G u G 2 , G 3 , G 4 )e lK ° x + c.c. ) 



(2tv/k 



,-i 



2tt/kq 



+ (2tt/k 



dx j dy { v* m d + q* m G 2 + ^* b0 \ z=0 G 3 } 

POO _ 

dy / dz ^* m G A , (3.42) 



"2tt/ko r x 

1 ~ 1 ' dx 



where C b0 = C bKo , v b0 = v bKo , qw = qb Ko , and \I>6o = ^b Ko - Completely expanding Eqn. |3.42 
gives the amplitude equation 



F 1 d T A + F 2 A + F 3 {2iK d x fA Q + F 4 A \A \ 2 = 



(3.43) 



in the slow scales X and T. The coefficients Fj are 



1/2 



Fi = / dy q* m q 

-1/2 
-1/2 



/lf 2 
dy{ -KlTZ c0 v* bo (qo - Q^oU=o) } 
-1/2 



(3.44) 
(3.45) 
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F 3 = f 1/2 dy{ {2K Q )- 1 v* w [2{d 2 - kI)% - n c0 % + n c0 Q^' \ z=0 } - v* bo v + (2K y 1 q* b0 ^' \ z=0 } 

J -1/2 

/OO /"OO _ _ 

dy / dz (2ko)- 1 ^^ ^ , (3.46) 
-OO J —OO 

Fa = f ^ dy{ (t Ko ) 2 n c0 v* bo [-q* (d y ^l\ z=0 ) - 2(d y q* )^{\z= + 2ql(d y %\ z=0 ) 
+ (d y ql)%\ z=0 + q {d y %\ z=0 ) - (d y q e 2 )y \z=o)} 

+Q*w[l(9 y vt)q* + (2iK )u* o ql + + v\{d y <&) + v (d y ql)] } , (3.47) 



where the prime denotes d K . In terms of the fast variables x and t, Eqn. |3.43| is expressed as 



r d t A = eA + - g A\A\ 2 , (3.48) 

such that A{x,t) = e l > 2 A (X,T), r = -F 1 /F 2 , Co = -4^F 3 /F 2 , and g = -F 4 /F 2 . 

The normalization of the amplitude function in the solution of the e 1//2 equation, 
Eqn. |3.10| , is arbitrary. The scale of A(x, t) can be set by requiring 



Nu - 1 = (qv)/(a s E y ) = \A\ 2 , (3.49) 

where 

i i-2-k/ko rl/2 

(...) = (271/ko)- 1 / dx / dy{...) . (3.50) 

JO J -1/2 

Note that Nu is the "Nusselt number" for the electroconvection problem, which is defined to 
be the ratio of the total current density to the conducted current density, spatially averaged. 

IV. DISCUSSION AND CONCLUSION 



To find the coefficients of the amplitude equation, Eqn. |3.48| , we evaluate Eqns. |3.44 - 



|3.47| using the numerical techniques described in Refs. and [|T7|]. The y-integrations are 
performed by the Romberg method. A simple SOR algorithm is employed to solve the 
Helmholtz equations, Eqns. |3.17| , |3.21| , and |B11| , on a N x N grid in the first quadrant 



of the yz plane. The solutions in the rest of yz plane follow from symmetry. The double 
integration in Eqn |3.46| is performed using a 2d trapezoidal rule based on the same grid. 
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The coefficients are extrapolated such that iV — > oo and Nfu m /N — > 0, where 2Nfu m is the 
number of grid points across the width of the film. The Fourier series in Eqns. [A4|, [A7|, 



and [B9] are expanded up to / = 29. Six modes are employed in the solutions Eqns. |2.25 - 



|2.27| , Eqn. and Eqns. |H5| - [FJ7| , of the linear stability, order e, and adjoint problems, 



respectively. Including more terms in these series expansions does not significantly change 
our final results. 

The values of To, £o> and go are shown in Table |. Only go, the coefficient of the nonlinear 
term, depends on the normalization of A by the Nusselt number according to Eqn. |3.49| . 
These results can be compared with those obtained independently from the linear stability 
calculations of Daya et al. || In the latter approach, the correlation length £ was derived 
from the curvature of the neutral curve at Ko and the characteristic time r from the linear 
growth rate at k . Both r and £o from the linear stability analysis are in good agreement 
with the present calculation. This provides a useful independent check of our numerical 
results. 

The comparison of our theoretical results with the experiments of Mao et al. [0] is difficult 
at the present time due to the uncertainties in the measurements of the material parameters 
of the liquid crystal film. For example, to nondimensionalize the experimentally measured 
value of go, the factor eo/cr 2 s 2 must be employed. While the thickness s of the smectic 
thin film can be measured accurately, the bulk conductivity of the liquid crystal is much 
less well characterized. Our calculated values of the (Nusselt normalized) value of go are 
substantially larger than those estimated from experiment , but in view of the uncertainty 
in u (roughly a factor of 3), no more precise comparison can currently be made. Experiments 
which will more accurately measure a and the viscosity rj in annular films are presently being 
performed. (T^j 

The flow pattern which develops just above onset can be visualized by evaluating the 
velocity vector field u on the xy plane. The lowest order x and y components of the 
field are given by Eqns. [A^ and |3.1(J| , respectively. The amplitude function A in these 
expressions is obtained by solving Eqn. |3.48| for the steady state case with a specified control 
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parameter e. An example of the resulting vortex pattern is shown in Fig. 0. This may 
be qualitatively compared with the experimental pattern shown in Fig. 6(b) of Ref. ||. 
As above, a quantitative comparison of theoretical and experimental velocities is difficult 
because of the experimental uncertainty in a. 

In conclusion, a multiple-scales expansion of the basic electrohydrodynamic equations 
for electroconvection in a suspended fluid film was used to find the lowest order amplitude 
equation. The set of basic equations were not self-adjoint, necessitating the evaluation of 
the adjoint eigenfunctions. The coefficients To, £o, and go of the resulting Ginzburg-Landau 
equation were determined by numerical integration. The results of this work can be employed 
in further studies of the weakly nonlinear phenomena near the onset of electroconvection in 
suspended smectic films. Of particular interest is the mechanism of wavelength selection |]J 
and effect of sidewalks on the convection pattern @,|14[ . 



ACKNOWLEDGMENTS 

We would like to thank John R. de Bruyn for numerous discussions. This research was 
supported by the Natural Sciences and Engineering Research Council of Canada. 

APPENDIX A: 



In this appendix, a method for solving the order e equations, Eqns. [3.14| - [3.2l| is described. 



Note that the functions v\(y) and are velocity fields and satisfy Eqn. |2.14j , the boundary 



conditions on v. The quantities q\(y) and q^iy) are charge densities. The functions ^\(y,z) 
and ty e 2 (y,z) are electric potentials and satisfy Eqns. |2.15| - |2.18| , the boundary conditions on 
*. 

To fully specify the order e equations, we relate the x- and ^-components of the velocity 
field. This is accomplished by expanding the incompressibility condition Eqn. p.2| via the 
multiple-scales method. At order e 1//2 , Eqn. pT2| is 
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d x uo(x, y) + d y v (x, y) = , (Al) 
which can be simplified using the order e 1 / 2 solution for v (x,y), Eqn. |3.10| , to 



u (x, y) = A (X, T)u (y)e tK0X + c.c. , (A2) 

where 

u (y) = in^idyVoiy)) . (A3) 



The solution of the first set of order e equations, Eqns. |3.14| - |3.17| , begins with Eqn. |3. 1 5 . 



In the linear stability calculation ||, the functions vo(y) and qo{y) are chosen to be even. 
So d y v (y) and d y q (y) are odd. The derivative of the base state charge density Q(y) is 
even which means that if v{(y) is odd, then the product Q{y)v\{y) is also odd. Hence the 
nonhomogeneous part of Eqn. |3.15| can be expanded in a Fourier sine series 



- (d y v (y))q (y) + v (y)(d y q (y)) + Q{y)v{{y) = J^ksin^irly) . (A4) 

i=i 

The general solution of Eqn. |3.15 



is 



%(y, 0) = Msinh{2K y) - ^[(2ttZ) 2 + (2KQ) 2 }-\sin{2irly) . (A5) 

i=i 

The coefficient M is zero due to the boundary condition ^(|,0) = 0. Assuming a trial 
solution 

N 

v{{y) = E m S m (y) , (A6) 

m=l 



where S m is the odd Chandrasekhar function [16 with E m = 1 for m = 1,...N, yields 



the coefficients b\ in Eqn. [A3]. Equation |A5| and the boundary conditions Eqns. [2.16| - [2.18| 



determine *ff\(y,z) via the Helmholtz equation, Eqn. |3.17| , which is solved by a numerical 
relaxation method. The function q\{y) is found by numerical differentiation of ^{(y, z) 
in Eqn. |3.16| . Then q\{y) and ty\(y,0) are substituted into Eqn. |3.14| to calculate a new 
estimate of v\(y). This process is repeated until v\{y), q{(y), and ^f[(y, z) are self-consistently 
determined. 



In the second set of order e equations, Eqns. |3.18| - |3.2l| , Eqn. |3.18| and the boundary 
conditions Eqn. [2.14| indicate that v\(y) = 0. The right-hand side of Eqn. |3.19| , simplified 
via Eqn. |A3| , is an odd function and can be expanded in a Fourier sine series. The general 
solution of Eqn. |3.19| is 



(A7) 



1=1 



The variable \P§(y, z) is specified by solving the Laplace equation, Eqn. p.21|, by a relaxation 
method, subject to Eqn. |A7| and the boundary conditions Eqns. [2.16| - [2.18[ The function 
q^(y) is numerically calculated via Eqn. |3.20| . 



APPENDIX B: 



The solution of the adjoint problem, Eqns. |3.28j - |3.37| , is discussed in this section. Sub- 
stitution of the vector C bKQ , Eqn. |3.38| , into Eqns. |3.29H3.3l2j gives 



(d y - k ) v bK0 - Qq bKo = 



-KllZ c0 Qv bK0 + (dy - nl)q bK0 + (d z ^ bK0 )\ z=0 + - {d z $ bK0 )\ z=0 - = 



(fl{ + ^-^)* teo =0 



(Bl) 
(B2) 
(B3) 
(B4) 



To simplify the notation, let v b0 = v bK0 , q b0 = q bK0 , and * 60 = * 6ft0 . The functions v b0 (y), 
q~w(y), and ^/ b o(y,z) are expanded as 

oo 

Vbo(y) = Yl BmVWmiy) , (B5) 
m=l 

oo 

?6o (y) = Yl B mqbom(y) , (B6) 

m=l 

oo 

f w (y, z) = B m ^ b0m (y, z) , (B7) 

m=l 

where t>&o m (2/)) q b om(y)-> an d ^b0m(y, z ) satisfy the adjoint boundary conditions, Eqns. |3.33 - 



ET371 
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With the solutions Eqns. |B5|-|B7|, Eqn. |B1| implies that 



(d y - k ) v Wm (y) - Q(y)q b0m (y) = 



(B8) 



Since q b o m {y) must satisfy q b0m {y = ±|) = 0, let qwm(y) = cos((2m - ljny). The product 
of Q(y) and cos ((2m — l)ivy) is even and can be represented by a Fourier cosine series. The 
general solution of Eqn. [B8| is 

oo 

VbOm(y) = Micosh(K y) + M 2 ysinh(n y) + ^[( 27r 2 + nl]~ 2 b m icos(2Trly) , (B9) 

1=0 

where the constants Mi and M 2 are specified by the boundary conditions v b o m (y = ±§) = 
(d y vwm){y = ±|) = to be 



M x = -2(«o + sm/i(/t ))" 1 [sm/i(K /2) + (k /2)cos/i(k /2)] ^(-1)' [(2ttZ) 2 + k 2 ]~ 2 & 



VrtZ 



M 2 = -2(«o + sinh(K ))- l [-K sinh(K Q /2)} ^(-l) i [(2 7 r/) 2 + k 2 ]- 2 & w 



Substitution of Eqns. |B5| - |B7| into Eqn. [B2| and Eqn. [B4] yields 



^bOmiV, 0) = -«o^-c0^bOm(2/) 



(BIO) 



and 



(d 2 y + d 2 z -K 2 )^ b0m (y,z) = o. 



(Bll) 



The latter is a Helmholtz equation subject to the adjoint boundary conditions Eqns. |3.36 
and |3.37| , with ^bOm(y, 0), \y\ < |, given by Eqn. |B10| . Equation |B11| is solved by a numerical 
relaxation method. 

Using Eqns. |B5|-|B7] to expand Eqn. IBB] leads to 



Bml-KlKcoQvwm + (dy - K 2 )q b0m + (d z ^ b0 m)\z=0+ ~ (d**60m)|*=0-] = . (B12) 



Since v b Q m (y), q b o m (y), and ^wm(y, z) are known functions, Eqn. |B12j implies that the coef- 
ficients B m vanish unless the compatibility condition 
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1/2 2 - - ~ 

dy[—K 'R-coQ(lbOnVbOm + QbOn(d ~ K o)(jbOrn + 2^0n (c^&Om) 1 2=0+] 
-1/2 y 



o, 

(B13) 



is satisfied. Note that the relation (d z ^? b om)\z=o- = ~ (9z^b0m)\ z =o+i which is analogous to 
the discontinuity in the electric field d z ^>o m across the two surfaces of the film in the direct 
problem, is used to derive Eqn. BT3|. In Eqn. IBTj , the values of k,q and TZ c q are fixed to be 
those obtained from the linear stability analysis. The coefficients B m are determined by the 
matrix equation 



nm-^m 



o 



(B14) 



where T nm is given in Eqn. |B13| . These coefficients are then substituted into Eqns. p5| - p7 
to generate Vbo(y), (jbo(y), and ^?bo(y, z )- This specifies Cb Ko and the general solution vector 
C b , Eqn. ggg . 
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TABLES 
TABLE I. Numerical results 





Nonlinear Analysis 


Linear Stability 


Wire Electrode Geometry 


critical wavenumber, k,q 


4.7467 


4.744 


critical control parameter, 7Z c q 


76.855 


76.77 


correlation length, £o 


0.28484 


0.2843 


time scale, To 


0.35072 


0.351 


nonlinear coupling, go 


1.74602 




Plate Electrode geometry 


critical wavenumber, no 


4.2239 


4.223 


critical control parameter, 7Z c o 


91.855 


91.84 


correlation length, £o 


0.29743 


0.2975 


time scale, to 


0.37155 


0.372 


nonlinear coupling, go 


2.8424 
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FIGURES 




•VI 2 



X 



FIG. 1. Schematic of film geometry and coordinate system. The wire electrode configuration 
is shown. The vortex pair periodicity is A = 2nd/ k, where d is the film width. The thickness s of 
the film (not shown) is such that s <C d. 
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FIG. 2. Vortex pattern just above onset. The dimensionless velocity must be scaled by sa/eo, 
where s, a, and eo are the film thickness, bulk conductivity, and permittivity of free space, respec- 
tively. Here we plot the vector velocity field for wire electrodes with control parameter e = 0.1. 
Using s = 142 nm and a = 2.0 x 10 - ' (fim) -1 , which are typical values for smectic films, gives 
sa/eo = 3.2 mm/s. The magnitude of the velocity at the centre of the figure is approximately 
equal to 2.4 mm/s. 
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